Rapid simulation of glycoprotein structures by grafting and steric exclusion of glycan conformer libraries.

Most membrane proteins are modified by covalent addition of complex sugars through N- and O-glycosylation. Unlike proteins, glycans do not typically adopt specific secondary structures and remain very mobile, shielding potentially large fractions of protein surface. High glycan conformational freedom hinders complete structural elucidation of glycoproteins. Computer simulations may be used to model glycosylated proteins but require hundreds of thousands of computing hours on supercomputers, thus limiting routine use. Here, we describe GlycoSHIELD, a reductionist method that can be implemented on personal computers to graft realistic ensembles of glycan conformers onto static protein structures in minutes. Using molecular dynamics simulation, small-angle X-ray scattering, cryoelectron microscopy, and mass spectrometry, we show that this open-access toolkit provides enhanced models of glycoprotein structures. Focusing on N-cadherin, human coronavirus spike proteins, and gamma-aminobutyric acid receptors, we show that GlycoSHIELD can shed light on the impact of glycans on the conformation and activity of complex glycoproteins.

IHMCIF: An Extension of the PDBx/mmCIF Data Standard for Integrative Structure Determination Methods.

IHMCIF (github.com/ihmwg/IHMCIF) is a data information framework that supports archiving and disseminating macromolecular structures determined by integrative or hybrid modeling (IHM), and making them Findable, Accessible, Interoperable, and Reusable (FAIR). IHMCIF is an extension of the Protein Data Bank Exchange/macromolecular Crystallographic Information Framework (PDBx/mmCIF) that serves as the framework for the Protein Data Bank (PDB) to archive experimentally determined atomic structures of biological macromolecules and their complexes with one another and small molecule ligands (e.g., enzyme cofactors and drugs). IHMCIF serves as the foundational data standard for the PDB-Dev prototype system, developed for archiving and disseminating integrative structures. It utilizes a flexible data representation to describe integrative structures that span multiple spatiotemporal scales and structural states with definitions for restraints from a variety of experimental methods contributing to integrative structural biology. The IHMCIF extension was created with the benefit of considerable community input and recommendations gathered by the Worldwide Protein Data Bank (wwPDB) Task Force for Integrative or Hybrid Methods (wwpdb.org/task/hybrid). Herein, we describe the development of IHMCIF to support evolving methodologies and ongoing advancements in integrative structural biology. Ultimately, IHMCIF will facilitate the unification of PDB-Dev data and tools with the PDB archive so that integrative structures can be archived and disseminated through PDB.

Structure and assembly of a bacterial gasdermin pore.

In response to pathogen infection, gasdermin (GSDM) proteins form membrane pores that induce a host cell death process called pyroptosis1-3. Studies of human and mouse GSDM pores have revealed the functions and architectures of assemblies comprising 24 to 33 protomers4-9, but the mechanism and evolutionary origin of membrane targeting and GSDM pore formation remain unknown. Here we determine a structure of a bacterial GSDM (bGSDM) pore and define a conserved mechanism of pore assembly. Engineering a panel of bGSDMs for site-specific proteolytic activation, we demonstrate that diverse bGSDMs form distinct pore sizes that range from smaller mammalian-like assemblies to exceptionally large pores containing more than 50 protomers. We determine a cryo-electron microscopy structure of a Vitiosangium bGSDM in an active 'slinky'-like oligomeric conformation and analyse bGSDM pores in a native lipid environment to create an atomic-level model of a full 52-mer bGSDM pore. Combining our structural analysis with molecular dynamics simulations and cellular assays, our results support a stepwise model of GSDM pore assembly and suggest that a covalently bound palmitoyl can leave a hydrophobic sheath and insert into the membrane before formation of the membrane-spanning ß-strand regions. These results reveal the diversity of GSDM pores found in nature and explain the function of an ancient post-translational modification in enabling programmed host cell death.

Structure and Internal Dynamics of Short RNA Duplexes Determined by a Combination of Pulsed EPR Methods and MD Simulations.

We used EPR spectroscopy to characterize the structure of RNA duplexes and their internal twist, stretch and bending motions. We prepared eight 20-base-pair-long RNA duplexes containing the rigid spin-label Çm, a cytidine analogue, at two positions and acquired orientation-selective PELDOR/DEER data. By using different frequency bands (X-, Q-, G-band), detailed information about the distance and orientation of the labels was obtained and provided insights into the global conformational dynamics of the RNA duplex. We used 19F Mims ENDOR experiments on three singly Çm- and singly fluorine-labeled RNA duplexes to determine the exact position of the Çm spin label in the helix. In a quantitative comparison to MD simulations of RNA with and without Çm spin labels, we found that state-of-the-art force fields with explicit parameterization of the spin label were able to describe the conformational ensemble present in our experiments. The MD simulations further confirmed that the Çm spin labels are excellent mimics of cytidine inducing only small local changes in the RNA structure. Çm spin labels are thus ideally suited for high-precision EPR experiments to probe the structure and, in conjunction with MD simulations, motions of RNA.

Mechanism of proton-powered c-ring rotation in a mitochondrial ATP synthase.

Proton-powered c-ring rotation in mitochondrial ATP synthase is crucial to convert the transmembrane protonmotive force into torque to drive the synthesis of adenosine triphosphate (ATP). Capitalizing on recent cryo-EM structures, we aim at a structural and energetic understanding of how functional directional rotation is achieved. We performed multi-microsecond atomistic simulations to determine the free energy profiles along the c-ring rotation angle before and after the arrival of a new proton. Our results reveal that rotation proceeds by dynamic sliding of the ring over the a-subunit surface, during which interactions with conserved polar residues stabilize distinct intermediates. Ordered water chains line up for a Grotthuss-type proton transfer in one of these intermediates. After proton transfer, a high barrier prevents backward rotation and an overall drop in free energy favors forward rotation, ensuring the directionality of c-ring rotation required for the thermodynamically disfavored ATP synthesis. The essential arginine of the a-subunit stabilizes the rotated configuration through a salt bridge with the c-ring. Overall, we describe a complete mechanism for the rotation step of the ATP synthase rotor, thereby illuminating a process critical to all life at atomic resolution.

Encoding prior knowledge in ensemble refinement.

The proper balancing of information from experiment and theory is a long-standing problem in the analysis of noisy and incomplete data. Viewed as a Pareto optimization problem, improved agreement with the experimental data comes at the expense of growing inconsistencies with the theoretical reference model. Here, we propose how to set the exchange rate a priori to properly balance this trade-off. We focus on gentle ensemble refinement, where the difference between the potential energy surfaces of the reference and refined models is small on a thermal scale. By relating the variance of this energy difference to the Kullback-Leibler divergence between the respective Boltzmann distributions, one can encode prior knowledge about energy uncertainties, i.e., force-field errors, in the exchange rate. The energy uncertainty is defined in the space of observables and depends on their type and number and on the thermodynamic state. We highlight the relation of gentle refinement to free energy perturbation theory. A balanced encoding of prior knowledge increases the quality and transparency of ensemble refinement. Our findings extend to non-Boltzmann distributions, where the uncertainty in energy becomes an uncertainty in information.

Mechanism and cellular function of direct membrane binding by the ESCRT and ERES-associated Ca2+-sensor ALG-2.

Apoptosis linked Gene-2 (ALG-2) is a multifunctional intracellular Ca2+ sensor and the archetypal member of the penta-EF hand protein family. ALG-2 functions in the repair of damage to both the plasma and lysosome membranes and in COPII-dependent budding at endoplasmic reticulum exit sites (ERES). In the presence of Ca2+, ALG-2 binds to ESCRT-I and ALIX in membrane repair and to SEC31A at ERES. ALG-2 also binds directly to acidic membranes in the presence of Ca2+ by a combination of electrostatic and hydrophobic interactions. By combining giant unilamellar vesicle-based experiments and molecular dynamics simulations, we show that charge-reversed mutants of ALG-2 at these locations disrupt membrane recruitment. ALG-2 membrane binding mutants have reduced or abrogated ERES localization in response to Thapsigargin-induced Ca2+ release but still localize to lysosomes following lysosomal Ca2+ release. In vitro reconstitution shows that the ALG-2 membrane-binding defect can be rescued by binding to ESCRT-I. These data thus reveal the nature of direct Ca2+-dependent membrane binding and its interplay with Ca2+-dependent protein binding in the cellular functions of ALG-2.

Hierarchical Assembly of Single-Stranded RNA.

Single-stranded RNA (ssRNA) plays a major role in the flow of genetic information-most notably, in the form of messenger RNA (mRNA)-and in the regulation of biological processes. The highly dynamic nature of chains of unpaired nucleobases challenges structural characterizations of ssRNA by experiments or molecular dynamics (MD) simulations alike. Here, we use hierarchical chain growth (HCG) to construct ensembles of ssRNA chains. HCG assembles the structures of protein and nucleic acid chains from fragment libraries created by MD simulations. Applied to homo- and heteropolymeric ssRNAs of different lengths, we find that HCG produces structural ensembles that overall are in good agreement with diverse experiments, including nuclear magnetic resonance (NMR), small-angle X-ray scattering (SAXS), and single-molecule Förster resonance energy transfer (FRET). The agreement can be further improved by ensemble refinement using Bayesian inference of ensembles (BioEn). HCG can also be used to assemble RNA structures that combine base-paired and base-unpaired regions, as illustrated for the 5' untranslated region (UTR) of SARS-CoV-2 RNA.

Three-step docking by WIPI2, ATG16L1, and ATG3 delivers LC3 to the phagophore.

The covalent attachment of ubiquitin-like LC3 proteins (microtubule-associated proteins 1A/1B light chain 3) prepares the autophagic membrane for cargo recruitment. We resolve key steps in LC3 lipidation by combining molecular dynamics simulations and experiments in vitro and in cellulo. We show how the E3-like ligaseautophagy-related 12 (ATG12)-ATG5-ATG16L1 in complex with the E2-like conjugase ATG3 docks LC3 onto the membrane in three steps by (i) the phosphatidylinositol 3-phosphate effector protein WD repeat domain phosphoinositide-interacting protein 2 (WIPI2), (ii) helix α2 of ATG16L1, and (iii) a membrane-interacting surface of ATG3. Phosphatidylethanolamine (PE) lipids concentrate in a region around the thioester bond between ATG3 and LC3, highlighting residues with a possible role in the catalytic transfer of LC3 to PE, including two conserved histidines. In a near-complete pathway from the initial membrane recruitment to the LC3 lipidation reaction, the three-step targeting of the ATG12-ATG5-ATG16L1 machinery establishes a high level of regulatory control.

Tracing the substrate translocation mechanism in P-glycoprotein.

P-glycoprotein (Pgp) is a prototypical ATP-binding cassette (ABC) transporter of great biological and clinical significance.Pgp confers cancer multidrug resistance and mediates the bioavailability and pharmacokinetics of many drugs (Juliano and Ling, 1976; Ueda et al., 1986; Sharom, 2011). Decades of structural and biochemical studies have provided insights into how Pgp binds diverse compounds (Loo and Clarke, 2000; Loo et al., 2009; Aller et al., 2009; Alam et al., 2019; Nosol et al., 2020; Chufan et al., 2015), but how they are translocated through the membrane has remained elusive. Here, we covalently attached a cyclic substrate to discrete sites of Pgp and determined multiple complex structures in inward- and outward-facing states by cryoEM. In conjunction with molecular dynamics simulations, our structures trace the substrate passage across the membrane and identify conformational changes in transmembrane helix 1 (TM1) as regulators of substrate transport. In mid-transport conformations, TM1 breaks at glycine 72. Mutation of this residue significantly impairs drug transport of Pgp in vivo, corroborating the importance of its regulatory role. Importantly, our data suggest that the cyclic substrate can exit Pgp without the requirement of a wide-open outward-facing conformation, diverting from the common efflux model for Pgp and other ABC exporters. The substrate transport mechanism of Pgp revealed here pinpoints critical targets for future drug discovery studies of this medically relevant system.

Neighbor List Artifacts in Molecular Dynamics Simulations.

Molecular dynamics (MD) simulations are widely used in biophysical research. To aid nonexpert users, most simulation packages provide default values for key input parameters. In MD simulations using the GROMACS package with default parameters, we found large membranes to deform under the action of a semi-isotropically coupled barostat. As the primary cause, we identified overly short outer cutoffs and infrequent neighbor list updates that resulted in missed nonbonded interactions. Small but systematic imbalances in the apparent pressure tensor then induce unphysical asymmetric box deformations that crumple the membrane. We also observed rapid oscillations in averages of the instantaneous pressure tensor components and traced these to the use of a dual pair list with dynamic pruning. We confirmed that similar effects are present in MD simulations of neat water in atomistic and coarse-grained representations. Whereas the slight pressure imbalances likely have minimal impact in most current atomistic MD simulations, we expect their impact to grow in studies of ever-larger systems with coarse-grained representation, in particular, in combination with anisotropic pressure coupling. We present measures to diagnose problems with missed interactions and guidelines for practitioners to avoid them, including estimates for appropriate values for the outer cutoff rl and the number of time steps nstlist between neighbor list updates.

Cell-free biosynthesis combined with deep learning accelerates de novo-development of antimicrobial peptides.

Bioactive peptides are key molecules in health and medicine. Deep learning holds a big promise for the discovery and design of bioactive peptides. Yet, suitable experimental approaches are required to validate candidates in high throughput and at low cost. Here, we established a cell-free protein synthesis (CFPS) pipeline for the rapid and inexpensive production of antimicrobial peptides (AMPs) directly from DNA templates. To validate our platform, we used deep learning to design thousands of AMPs de novo. Using computational methods, we prioritized 500 candidates that we produced and screened with our CFPS pipeline. We identified 30 functional AMPs, which we characterized further through molecular dynamics simulations, antimicrobial activity and toxicity. Notably, six de novo-AMPs feature broad-spectrum activity against multidrug-resistant pathogens and do not develop bacterial resistance. Our work demonstrates the potential of CFPS for high throughput and low-cost production and testing of bioactive peptides within less than 24 h.

Mechanism and cellular function of direct membrane binding by the ESCRT and ERES-associated Ca2+-sensor ALG-2.

Apoptosis Linked Gene-2 (ALG-2) is a multifunctional intracellular Ca2+ sensor and the archetypal member of the penta-EF hand protein family. ALG-2 functions in the repair of damage to both the plasma and lysosome membranes and in COPII-dependent budding at endoplasmic reticulum exit sites (ERES). In the presence of Ca2+, ALG-2 binds to ESCRT-I and ALIX in membrane repair and to SEC31A at ERES. ALG-2 also binds directly to acidic membranes in the presence of Ca2+ by a combination of electrostatic and hydrophobic interactions. By combining GUV-based experiments and molecular dynamics simulations, we show that charge-reversed mutants of ALG-2 at these locations disrupt membrane recruitment. ALG-2 membrane binding mutants have reduced or abrogated ERES localization in response to Thapsigargin-induced Ca2+ release but still localize to lysosomes following lysosomal Ca2+ release. In vitro reconstitution shows that the ALG-2 membrane-binding defect can be rescued by binding to ESCRT-I. These data thus reveal the nature of direct Ca2+-dependent membrane binding and its interplay with Ca2+-dependent protein binding in the cellular functions of ALG-2.

Efficient generation of random rotation matrices in four dimensions.

Four-dimensional (4D) rotations have applications in the fields of robotics, computer vision, and rigid-body mechanics. In the latter, they can be used to transform between equimomental systems of point masses. Here we provide an efficient algorithm to generate random 4D rotation matrices covering an arbitrary, predefined range of rotation angles. These matrices can be combined with Monte Carlo methods for the efficient sampling of the SO(4) group of 4D rotations. The matrices are unbiased and constructed such that repeated rotations result in uniform sampling over SO(4). The algorithm can be used to optimize the mass partitioning in coarse-grained simulation models of molecules involving coupled constraints for stable time integration.

Martini 3 Coarse-Grained Force Field for Cholesterol.

Cholesterol plays a crucial role in biomembranes by regulating various properties, such as fluidity, rigidity, permeability, and organization of lipid bilayers. The latest version of the Martini model, Martini 3, offers significant improvements in interaction balance, molecular packing, and inclusion of new bead types and sizes. However, the release of the new model resulted in the need to reparameterize many core molecules, including cholesterol. Here, we describe the development and validation of a Martini 3 cholesterol model, addressing issues related to its bonded setup, shape, volume, and hydrophobicity. The proposed model mitigates some limitations of its Martini 2 predecessor while maintaining or improving the overall behavior.

Molecular mechanisms of inorganic-phosphate release from the core and barbed end of actin filaments.

The release of inorganic phosphate (Pi) from actin filaments constitutes a key step in their regulated turnover, which is fundamental to many cellular functions. The mechanisms underlying Pi release from the core and barbed end of actin filaments remain unclear. Here, using human and bovine actin isoforms, we combine cryo-EM with molecular-dynamics simulations and in vitro reconstitution to demonstrate how actin releases Pi through a 'molecular backdoor'. While constantly open at the barbed end, the backdoor is predominantly closed in filament-core subunits and opens only transiently through concerted amino acid rearrangements. This explains why Pi escapes rapidly from the filament end but slowly from internal subunits. In a nemaline-myopathy-associated actin variant, the backdoor is predominantly open in filament-core subunits, resulting in accelerated Pi release and filaments with drastically shortened ADP-Pi caps. Our results provide the molecular basis for Pi release from actin and exemplify how a disease-linked mutation distorts the nucleotide-state distribution and atomic structure of the filament.

Structure-function analysis suggests that the photoreceptor LITE-1 is a light-activated ion channel.

Sensation of light is essential for all organisms. The eye-less nematode Caenorhabditis elegans detects UV and blue light to evoke escape behavior. The photosensor LITE-1 absorbs UV photons with an unusually high extinction coefficient, involving essential tryptophans. Here, we modeled the structure and dynamics of LITE-1 using AlphaFold2-multimer and molecular dynamics (MD) simulations and performed mutational and behavioral assays in C. elegans to characterize its function. LITE-1 resembles olfactory and gustatory receptors from insects, recently shown to be tetrameric ion channels. We identified residues required for channel gating, light absorption, and mechanisms of photo-oxidation, involving a likely binding site for the peroxiredoxin PRDX-2. Furthermore, we identified the binding pocket for a putative chromophore. Several residues lining this pocket have previously been established as essential for LITE-1 function. A newly identified critical cysteine pointing into the pocket represents a likely chromophore attachment site. We derived a model for how photon absorption, via a network of tryptophans and other aromatic amino acids, induces an excited state that is transferred to the chromophore. This evokes conformational changes in the protein, possibly leading to a state receptive to oxidation of cysteines and, jointly, to channel gating. Electrophysiological data support the idea that LITE-1 is a photon and H2O2-coincidence detector. Other proteins with similarity to LITE-1, specifically C. elegans GUR-3, likely use a similar mechanism for photon detection. Thus, a common protein fold and assembly, used for chemoreception in insects, possibly by binding of a particular compound, may have evolved into a light-activated ion channel.

Crosstalk between regulatory elements in disordered TRPV4 N-terminus modulates lipid-dependent channel activity.

Intrinsically disordered regions (IDRs) are essential for membrane receptor regulation but often remain unresolved in structural studies. TRPV4, a member of the TRP vanilloid channel family involved in thermo- and osmosensation, has a large N-terminal IDR of approximately 150 amino acids. With an integrated structural biology approach, we analyze the structural ensemble of the TRPV4 IDR and the network of antagonistic regulatory elements it encodes. These modulate channel activity in a hierarchical lipid-dependent manner through transient long-range interactions. A highly conserved autoinhibitory patch acts as a master regulator by competing with PIP2 binding to attenuate channel activity. Molecular dynamics simulations show that loss of the interaction between the PIP2-binding site and the membrane reduces the force exerted by the IDR on the structured core of TRPV4. This work demonstrates that IDR structural dynamics are coupled to TRPV4 activity and highlights the importance of IDRs for TRP channel function and regulation.

Structure and assembly of a bacterial gasdermin pore.

In response to pathogen infection, gasdermin (GSDM) proteins form membrane pores that induce a host cell death process called pyroptosis1-33. Studies of human and mouse GSDM pores reveal the functions and architectures of 24-33 protomers assemblies4-9, but the mechanism and evolutionary origin of membrane targeting and GSDM pore formation remain unknown. Here we determine a structure of a bacterial GSDM (bGSDM) pore and define a conserved mechanism of pore assembly. Engineering a panel of bGSDMs for site-specific proteolytic activation, we demonstrate that diverse bGSDMs form distinct pore sizes that range from smaller mammalian-like assemblies to exceptionally large pores containing >50 protomers. We determine a 3.3 Å cryo-EM structure of a Vitiosangium bGSDM in an active slinky-like oligomeric conformation and analyze bGSDM pores in a native lipid environment to create an atomic-level model of a full 52-mer bGSDM pore. Combining our structural analysis with molecular dynamics simulations and cellular assays, our results support a stepwise model of GSDM pore assembly and suggest that a covalently bound palmitoyl can leave a hydrophobic sheath and insert into the membrane before formation of the membrane-spanning ß-strand regions. These results reveal the diversity of GSDM pores found in nature and explain the function of an ancient post-translational modification in enabling programmed host cell death.

Unwrapping NPT Simulations to Calculate Diffusion Coefficients.

In molecular dynamics simulations in the NPT ensemble at constant pressure, the size and shape of the periodic simulation box fluctuate with time. For particle images far from the origin, the rescaling of the box by the barostat results in unbounded position displacements. Special care is thus required when a particle trajectory is unwrapped from a projection into the central box under periodic boundary conditions to a trajectory in full three-dimensional space, e.g., for the calculation of translational diffusion coefficients. Here, we review and compare different schemes in use for trajectory unwrapping. We also specify the corresponding rewrapping schemes to put an unwrapped trajectory back into the central box. On this basis, we then identify a scheme for the calculation of diffusion coefficients from NPT simulations, which is a primary application of trajectory unwrapping. In this scheme, the wrapped and unwrapped trajectory are mutually consistent and their statistical properties are preserved. We conclude with advice on best practice for the consistent unwrapping of constant-pressure simulation trajectories and the calculation of accurate translational diffusion coefficients.